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Abstract 

We describe how an approximation to the minimal energy Skyrme 
crystal can be obtained from the holonomy of a Yang-Mills instanton. 
The appropriate instanton is twisted on a four-torus and has instanton 
number equal to one half. It generates a Skyrme field with the correct 
topological and symmetry properties of the crystal. An explicit solu- 
tion for the instanton is not known, but an analytical fit to numerical 
data is available and using this we obtain a Skyrme crystal whose 
energy is only 2% above that of the (numerically) known solution. 
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1 Introduction 



The Skyrme model is an effective theory of low energy hadron physics in 
which nucleons are described by topological solitons. The baryon number 
B is identified with the topological charge (or number of solitons) and it is 
of interest to determine the classical soliton solutions, since they describe 
the classical state of the nucleons. However, the Skyrme field equation is a 
nonlinear partial differential equation which is difficult to solve. In fact, even 
for the simplest case of B = 1 the minimal energy solution is only known 
numerically. 

In a region with high baryon number, the solitons form a crystal. Kle- 
banov |IJ began the study of Skyrme crystals and over the last 10 years 
various numerical and analytical techniques have been applied to find the 
solution with minimum energy per baryon. 

It has proved useful for understanding the solitons of the Skyrme model 
to construct approximate solutions by calculating the holonomy of Yang- 
Mills instantons ||. This method has led to a good understanding of the 
geometry of the B = 2 sector [^J, and the tetrahedral and cubic symmetries 
of the minimal energy B = 3 and B = 4 solutions respectively ||. The 
instanton method has computational advantages, and it provides a means by 
which the full soliton theory may be truncated to a finite number of degrees 
of freedom. However, one area in which it has yet to make a contribution is 
that of the Skyrme crystal. In this letter we take the first steps in applying 
the instanton method to the Skyrme crystal by identifying a suitable Yang- 
Mills instanton. We also make a qualitative calculation based on numerical 
data for the instanton. 

2 Skyrme Crystal 

The field U(x) of the static Skyrme model is an SU(2) matrix-valued function 
of the space coordinates x = (x\,X2,xs) = (x,y,z). It may be expressed in 
terms of sigma and pion fields as 

U = a + iir • r (1) 

where 7r = (tti, tt 2 , vr 3 ) is a triplet of pion fields and r denotes the Pauli 
matrices. The normalization constraint is a 2 + 7r.7r = 1. For finite energy 
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fields, it is required that as |x| — ► oo, C/(x) — > 1. The energy density is 

£ = -^~2 T < R ^) ~ ^TrdR^R^R,,^}) (2) 
and the baryon density is 

B = -^e ijk Tr{R i R j R k ) (3) 

where Ri is the right current dJJU^ 1 . The spatial integrals of £ and B give 
the energy E and the baryon number (topological charge) B of the Skyrme 
field. The above units are chosen so that the Fadeev-Bogomolny bound on 
the energy E is simply 

E > \B\. (4) 

In flat space this bound cannot be attained (except by the vacuum solution 
U = 1) and the single Skyrmion (B = 1) solution has energy E = 1.23. 

With certain relative orientations well separated Skyrmions attract, and 
at high density form a crystal. Klebanov's original crystal [|IJ was a simple 
cubic lattice of Skyrmions whose symmetries were motivated by the attempt 
to rotate the Skyrmions so as to maximize the attraction between nearest 
neighbours. It has energy per baryon of E = 1.08. Other symmetries were 
proposed which lead to slightly lower, but not minimal energy crystals ||. 
Following the work of Castillejo et al. ||, and of Kugler and Shtrikman [J/J, it 
is now understood that it is best to arrange the Skyrmions initially as a face- 
centred cubic lattice. Their orientations can be chosen very symmetrically to 
give maximal attraction between all nearest neighbours. This configuration 
relaxes to the minimal energy crystal, the Skyrme crystal, and in the process 
the symmetry increases further. The numerical investigations of , and the 
variational approach of using fourier series, lead to the conclusion that 
the Skyrme crystal has energy per baryon of E = 1.038 and a cubic unit cell 
of side L — 4.7. It is a crystal of half-Skyrmions, as a unit cell contains half 
a unit of baryon number. The fields are strictly periodic after translation by 
2L in the x,y or z directions, and within a cube of side 2L the baryon number 
is B = 4. The crystal symmetries, which differ from those of the crystals 
proposed in refs. [1] and [5], involve translations, reflections and rotations of 
the space coordinates combined with 0(4) rotations of the sigma and pion 



3 



fields. Explicitly, the symmetry generators are J7J 



(5) 
(6) 
(7) 



The fields obtained numerically, or by optimising the fourier series, are 
very well approximated by the analytical formulae 

a = cic 2 c 3 
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where 

Sj = sin(— -) and c$ = cos(— -A). (10) 



3 Twisted Instantons on T 4 

The pioneering work on Yang-Mills instantons on a four-torus T 4 with twisted 
boundary conditions was done by 't Hooft ||, motivated by the study of 
quark confinement in QCD. Let be the su(2)-valued Yang-Mills gauge 
potential, with field strength F^ v . Here greek indices take the values to 3, 
and latin indices 1 to 3, and we write x^ = (t,x,y,z), with t the euclidean 
time coordinate. We model the four-torus by the euclidean box defined by 
< x^ < a M . To give the boundary conditions we introduce Si7(2)-valued 
twist matrices f2 M each of which is independent of the /ith spacetime coor- 
dinate. The twisted boundary conditions are that the gauge potentials are 
periodic modulo gauge transformations 

A v {x ll = o M ) = n- x A v { Xll = 0)^ + n- x d v %. (11) 

The notation A u (x^ = a M ) means that the /ith coordinate is set equal to 
but the remaining coordinates are arbitrary. On 2-faces of the box the 
compatibility conditions which arise from QTTJ) are 

tt^Xv = a^Q^Xf, = a M ) = Vt v (x^ = a^)Vt^(x u = a u )Z^ u (12) 
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where each Z^ v is an element of the centre 7L.2 of SU(2). We introduce the 
gauge invariant antisymmetric twist integers n Mi , G Z (mod 2) by 

Zpv = exp(i7rn MI/ ). (13) 

It is these that are the important quantities, since any twist matrices with 
the same set of twist integers are gauge equivalent. It is convenient to collect 
the twist integers into a magnetic flux vector m and an electric flux vector 
k, where riij = e^ktrik and fc, = n Qi . 

't Hooft pointed out that the usual expression for the instanton number 
does not have to be integer valued in the presence of twist. In fact 

N =Y^L ^F^n** = 9 + ^ (14) 

where q is an integer. The mathematical explanation of this result can be 
found in |J , [[l(J . The expression ( |14|) is minus the second Chern number for 
an SU (2) gauge theory, but since the gauge functions are only defined modulo 
the centre of the gauge group, we are in fact dealing with the gauge group 
SU{2)/7L2 = SO(3), and that is why iV can be an integer or half-integer. 

Many elegant mathematical results are known for twisted instantons on 
T 4 , such as an existence proof Jll), a (non-explicit) construction of the appro- 
priate bundles in terms of K-theory [IIJ and a Nahm duality transformation 



The moduli (parameter) space of instantons of instanton number N is 



8N, that is, an integer multiple of four fTB"]. 



4 Skyrme Crystal from an Instanton 

We now show how an approximation to the Skyrme crystal described in 
section 2 can be obtained from the holonomy of a suitable instanton on T 4 . 

In order to obtain a Skyrme field with the correct cubic symmetry we 
set di = L and a = T, and choose isotropic magnetic and electric flux 
vectors m = (1,1,1) and k = (1,1,1). Then (|H|) becomes N = q + §. 
There are explicit abelian solutions || for iV = |, when the four-torus has 
the same length in each direction, but we are interested in approximating a 
crystal composed of half-Skyrmions so we wish to consider the case N — ^. 
Note that the Yang- Mills action on T 4 is bounded by 87r 2 |iV|, so the N — ^ 
instantons are the minimal action configurations outside the vacuum sector. 
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Moreover, their moduli space is four-dimensional and simply parametrises the 
translates, and certain discrete gauge transformations of, a unique N — \ 
instanton on the given torus. This instanton is invariant under the cubic 
subgroup of the spatial rotation group, and may be assumed to be centred 
with its maximum action density aix = y = z = L/2 and t = T/2. 

By choosing explicit twist matrices we perform a partial gauge fixing. 
Here it is convenient to use the Pauli matrices for the spatial twist matrices 

fij = ITj. (15) 

It is a simple matter to check that this choice satisfies the spatial components 
of (|12D with the twist integers given earlier. It is also easy to construct a 
twist matrix Q which satisfies the remaining components of (|T2|). Although 
we shall not need an explicit form for this matrix we give a possible choice 
for completeness, 

= (n a) + *v^5E^c.t (i6) 
t=i v c • c 

where q = cos(irxi/ L). 



The Skyrme field generated by the instanton is, by definition, the holon- 
omy of the instanton along all lines parallel to the euclidean time axis, ie 

C/(x) = (vexp £ A (x ,x) dx ^j fi„(x) (17) 

where V denotes path ordering. Note the inclusion of the twist matrix since 
the path must be a closed loop. The Skyrme field is defined on a cube of 
side L, and is extended by symmetry and continuity to a crystal with this 
cube as unit cell. The usual argument for instantons on 1R 4 can be repeated 
in this case to show that the instanton number N equals the baryon number 
B of the Skyrme field that it generates. Hence the Skyrme field has B — | 
in a unit cell. 

We can show, without using the explicit form (|T6|), that the Skyrme field 
has the symmetry (§) of a half-Skyrmion crystal. Using ( [TTf) and fll5|) we 
have 

U(x + L, y, z) = n (v exp A (x , x,y, z) dx^j n Vt Q (x + L,y, z). (18) 
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Equation (|T2| ) implies that 



Q (x + L,y,z) = -TiQ,o(x,y,z)Ti, (19) 

and combining (]I8|) and ( ]TD| ) gives 

U(x + L,y,z) = -TiU(x,y,z)n, (20) 

which is equivalent to symmetry (|8|). The symmetries (H)-(0) are a conse- 
quence of the cubic symmetry of the instanton. 

To progress, we would like to have an explicit analytical solution for the 
instanton of interest, but this is not available. Fortunately some numerical 
studies have been performed by Garcia Perez et al. |L4] and we can make 



use of their results as a first step in determining the instanton-generated 
Skyrme crystal, and estimating how good an approximation it is. Garcia 
Perez et al. treat the Yang-Mills theory as a lattice gauge theory and search 
for the minimal action configuration using the lattice cooling method. For 
the case T = oo they provide an analytical fit to the numerical data for the 
holonomy, and it is this result which we shall use here. Although T = oo 
the instanton is still localized in the euclidean time direction, with a scale 
which is determined by the spatial length of the four-torus, L. (In this limit, 
one lets the range of t be — oo < £ < oo and assumes that the instanton has 
its maximal action density at t — 0.) The analytical fit from jTJJ] for the 
holonomy in the euclidean time direction is 



U(x) = (II + i\l x " ( ^= A)2 f-r (21) 



where 



/ ii = cos(^)(l + asin 2 (^)), (22) 

/i = -sin(^)(l+/W(^)). (23) 

Here a = —0.196 and (5 = —0.172 are constants whose values are determined 
by the fit. Note that this Skyrme field has a similar form to the approximation 
to the Skyrme crystal (§) mentioned earlier. It is a simple matter to check 
explicitly that it has all the symmetry properties @ to (§|) of the Skyrme 
crystal. 
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It is relatively easy to obtain highly accurate numerical integrals for quan- 
tities such as the baryon number and energy of the Skyrme field (^l|) since 
the differentiation of U can be performed analytically, and furthermore the 
region of integration is finite. Integrating the baryon density (|3]) over a cube 
of side L confirms the baryon number to be B = |, and we find the energy 
per baryon is 

2 6146 

£ = 0.1069L + — — . (24) 

Lj 

Hence the minimum energy is E = 1.058 at L — 4.9, with both these figures 
being reasonably close to the Skyrme crystal values. One may take the view 
that the crystal size should be fixed at the value L = 4.7. The approximation 
( |2l"l) then contains no free parameters at the expense of an energy increase of 
0.1%. The fact that the energy is only 2% above that of the Skyrme crystal 
is encouraging, given that we are using only a fit to numerical data. Errors 
will be introduced both in the numerical construction of the instanton and 
in the analytical fitting process, and it would be very interesting if a more 
accurate numerical method could be employed in calculating the instanton 
and its holonomy. 

In [|TJ|] it is noted that moderate variations of the parameters a and j3 
may not significantly alter the quality of the fit to the numerical data. To see 
whether a more accurate scheme may produce a better approximate Skyrme 
crystal, we allow a and (3 to be arbitrary parameters. Using a steepest decent 
algorithm in (a, (3) space we find that the minimum crystal energy is obtained 
when a = —0.02 and f3 = —0.25. At these values the energy is E = 1.040 for 
a crystal size L = 4.7, ie the energy is within |% of the true minimum. 

Allowing T to vary could produce a better approximation to the Skyrme 
crystal. In [[Tj| it was shown that a slightly better approximation to the 
B = 1 Skyrmion in R 3 can be obtained from an instanton on 1R 3 times a finite 
euclidean time interval (a caloron) than from an instanton on R 4 . Moreover, 
a lower dimensional analogue of the procedure described in this letter exists 
[|16|j , where sine-Gordon kink chains are approximated by the holonomies 
of (DP 1 instantons on T 2 . In that example a much better approximation 
is obtained by having a finite length for the torus in the euclidean time 
direction. If this analogue result is a good guide to what happens in the case 
of the Skyrme crystal then one must certainly consider finite T. 
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5 Conclusion 



In this letter we have identified the Yang-Mills instanton whose holonomy 
produces a Skyrme field which has the topological and symmetry properties 
of the Skyrme crystal. Since an explicit solution is not known for the in- 
stanton it is not an easy task to obtain accurate quantitative information 
about the generated crystal. However, by using numerical work of Garcia 



Perez et al. [14] we have obtained some initial results which are encouraging. 



We hope that a more detailed investigation of this problem will prove fruitful. 
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